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ABSTRACT 


Tests  for  exponential  versus  IFRA  distributions 
based  on  incomplete  data  are  defined  and  shown  to 
be  unbiased.  The  tests  are  motivated  by  a  class 
of  tests  considered  in  detail  by  Bickel  and  Doksum. 
Tests  for  exponential  versus  IFR  distributions 
based  on  the  ranks  of  total  time  on  test  statistics 
are  also  considered. 


A  NOTE  ON  TESTS  FOR  MONOTONE 
FAILURE  RATE  BASED  ON  INCOMPLETE  DATA 

by 

Richard  E.  Barlow^  and  Frank  Proschan 


1.  INTRODUCTION  AND  SUMMARY 

Let  0  =  X/ri.  <  X.,  s  <  X,-.  <  ...  <  X,  .  be  the  order  statistics  of  a 
(0)  -  (1)  -  (2)  -  -  (n) 

(complete)  random  sample  from  a  population  with  distribution  F  and  density  f 
such  that  F(0)  *  0  .  Bickel  and  Doksum  (1968)  consider  the  problem  of  testing 


H, 


F(t)  -  1  -  e 


-At 


t  >  0  ,  A  >  0 


versus 


Hx  :  F  IFR 


(i.e.,  -log[l  -  F ( t ) ]  convex  on  [0,®))  . 

Let  D^  -  (n-i+l)(X^  -  X^  ^ )  ,  i  -  1,2,  ...»  n  .  They  consider  tests 
based  on  statistics  of  the  form 


l  aiDi 

i-1 

n 

l  »! 

i-1 


when 


where  a.  >  a„>  . . .  >  a  .  The  test,  $  ,  rejects  Hn 

1—2—  —  n  a  J  0 

n  /  n 

J  a  D  /  [  D  >  c  .  They  compute  the  asymptotic  relative  efficiency  of 

i-i  1  y  i-i  1  a,a,n 
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various  such  tests  relative  to  selected  parametric  alternatives.  Such  tests  were 
shown  to  be  unbiased  against  IFRA  (for  increasing  failure  rate  average) 
alternatives  by  Barlow  and  Proschan  (1966)  and  hence  a  fortiori  for  IFR 
alternatives.  [See  also  Birnbaum,  Esary  and  Marshall  (1965)  for  justification  of 
the  IFRA  assumption.) 

The  purpose  of  this  note  is  to  show  that  analogous  tests  designed  to  treat 
incomplete  samples  of  failure  data  are  also  unbiased  against  IFRA  alternatives. 
Let  be  the  time  to  failure  of  the  iC^  item  in  a  sample  of  size  n  .  Let  L^ 

be  a  given  truncation  time  for  the  i*”^  item  and  let 

*  min(Xi»Li)  i  *  1,2,  ....  n  . 


Let  0  =  Z/r..  <  Z,,.  <  ...  <  Z.,  .  be  the  first  k  observed  failure  times.  Note 
that  "withdrawals"  may  occur  between  Z^  and  Z(i+i)  anc*  t^iat  ^  is,  *n 
general,  a  random  variable.  Let  n(u)  be  the  (random)  number  of  items  on  test  at 
time  u  . 

* 

We  define  a  test,  4>a  ,  (a  modification  of  $a)  which  rejects  in  favor  of 

i 

:  F  IFRA 


(i.e.,  -{log[l  -  F(t) 1 }/t  nondecreasing  on  (O,®))  when 


r 


st 


Note  that  /  n(u)du  represents  the  total  time  on  test  between  the  i-1  and 

(i-l) 
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it*1  observed  failures.  The  distribution  of  can  be  computed  under  using 


the  fact  that  Y 


“(o’ 


n(u)du  (i  ■=  1,2,  ...»  k)  are  distributed  as 


( i— 1 ) 


independent  exponential  random  variables  under  Hq  conditioned  on  the  value  of  k 
* 

We  show  that  4>  is  an  unbiased  test  for  IFRA  alternatives  for  weights 

fl 

a  -  (a.,a„,  . a  )  for  which  a,  >  a_>  . . .  >  a 
i  i  n  l  —  i.  —  —  n 
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2.  DISTRIBUTION  OF  W  UNDER  H. 

_ a _ 0 

Let  r(t)  =  f(t)/[l  -  F ( t ) ]  be  the  failure  rate  function  for  F  We  will 
need  the  following  lemma,  stated  without  proof  in  Bray,  Crawford,  and  Proschan 
(1967). 

Lemma  1. 

For  any  distribution  F  (F (0 )  ■  0)  with  failure  rate  r(t)  , 

Z 


■.  •  r 

Z(i-1) 


r(u)n(u)du  ,  i  ■  1,2,  ...,  k  are  independently  distributed  with 


density  e  y  . 


Proof : 


Let  Y 


>(1)  r 

1  -  /  r(u)n(u)du  and  SqU)  “  /  r(u)n(u)du  .  Note  that  S^(t) 

0  0 

is  well  defined  up  to  the  time  of  the  first  observed  failure  since  n(u)  depends 
only  on  the  specified  truncation  times  L^  (i  *■  1,2,  . . . ,  n)  up  until  Zi 
Then 


(1)  * 


1  ’V  ‘  ‘“0 '*'(!) 7  '  yl]  "  P[L(1) 


P(Y,  >  yj  -  P[Sn(ZM.)  >  y. 


Z  / 1  v  ^  s 


exp[_so(so1<yi))]  ■  e  1 


Thus  Y^  has  density  e 

r(2)  r 

Now  let  ^2  "  /  r(u)n(u)du  and  (t)  *  J  r(u)n(u)du  .  Note  that 

Z(l)  1  X1 

conditionally  on  Z^  -  x^  ,  S*  is  well  defined  for  <_  t  <  Z ^  •  Hence 
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P[V2  “  y2  I  Z(l)  ‘  ’‘ll  *  P 


[Sx1(Z<2))  ' 


'2  1  "(1) 


P  Z(2)  >  ^Xj^2>  I  2 


(1)  '  “l]  '  eXp[-S*i(:\<y2))] 


-y- 


Thus  Y0  is  independent  of  and  also  exponentially  distributed  with  mean  1. 
If  we  continue  in  this  manner,  conditioning  on  previous  events,  we  establish  the 
lemma.  { | 

Under  ,  r(t)  =  X  and  we  see  from  the  lemma  that,  given  k  observed 
failures. 


W 

a 


k 

I 

st  i“l 


3iYi 


k 

l 

j-1 


I 


S  t 

where  **  denotes  stochastic  equality  and  Y^.Yj . Yk 

exponentially  distributed  random  variables  with  unit  mean. 


are  independent, 


♦ 


I 


3.  UNBIASEDNESS  UNDER  I  FRA  ALTERNATIVES 


We  need  the  following  lemma  to  establish  unbiasedness.  Define 


J *  r(u)du 


and  T(t) 


/" 


(u)du  . 


Lemma  2 . 


R ( t )  T (t) 

If  — — L  is  nondecreasing  in  t  ^  0  ,  n(t)  0  ,  and  • 11  “•  ^  is 

nonincreasing  in  t  >_  0  ,  then 


r  (u)d 


L 

u  f  r(u 


)dT(u) 


(i)  r(t)  > 


C 

J ”  r(u)dT(u) 


is  nondecreasing  in  t  >_  0  , 


when  the  indicated  integrals  exist. 


Proof : 


To  show  (i).  The  first  inequality  follows  from  differentiating 
R(t) 

Since  ~  ^  ^0  is  nondecreasing  in  t  >_  0  ,  we  can  approximate  R(t) 
arbitrarily  closely  from  below  by  a  positive  linear  combination  of  functions  of 
the  form 


0  0  <  t  <  x 


t  t  >  x 
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[cf.  Barlow,  Me^shall,  and  Proschan  (1967)].  By  the  Lebesque  monotone 
convergence  theorem,  we  need  only  establish  the  second  inequality  in  (ii)  for 
functions  R(t)  of  this  type.  Hence  for  t  >_  x  , 


v- 

J*  n(u)dR(u)  n(x)x  +  n(u)du 


This  is  nondecreasing  in  t  ^  x  since  (a)  T(t)  is  nondecreasing  in  t  0  , 

T  ( x ) 

and  (b)  xn(x)  -  T(x)  <_  0  since  is  nonincreasing  in  x  0  . 

To  show  (ii).  Clearly 


if  and  only  if 


(t)n(t)  J'  n(u)du  n(t)  J"  r(u)n(u)du 


which  follows  from  (i). 

Note  that  if  r(t)  is  nondecreasing  in  t  ^  0  ,  then  (ii)  follows  for  all 
n(t)  2.  0  • 

Lemma  2  may  be  used  in  testing  for  IFRA  in  models  other  than  the  one  described 
in  the  introduction;  see  for  example  the  model  of  Bray,  Crawford,  and  Proschan 
(1967). 


V 
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Theorem  1. 

If  F  is  IFRA  with  failure  rate  r(t)  and  Z...  <  •  ...  •  7..  .  are 

(1 )  —  ( .  )  —  -  (k) 

the  observed  failure  times,  n(t)  0  for  t  >_  0  ,  and  0  is  non  Increa  sing 

in  t  0  ,  then  (conditional  on  k) , 


k  t,U) 

l  -i  / 

i-1  1  J 


n(u)du 


liiziL 


/k)„ 


(u)du 


k 

,  A  *iYi 

at  1-1 _ 

k 

l  \ 

i-1  1 


where  a^  _>  i  •••  i  an  and  Y1,Y2’  ****  ^k  are  *IK*ePen<*ently  distributed  as 
exponential  random  variables  with  unit  mean. 

Proof: 

Since  n(u)  _>  0  and  T(t)/t  is  non increasing.  Lemma  2  applies,  yielding 


r 


i) 


r (u)n(u)du 


i  0 


a.  Z 


i  Ji) 

JL  " 


(u)du 


nondecreasing  in  i  -  1,2,  ...»  k  .  By  Lemma  1  we  need  only  show  that 


k  *(i) 

l  -i  / 

1-1 


n(u)du 


k 

I  ai  / 

i-1 


i) 


r (u)n(u)du 


(1) 


likJLL 


(l-t) 


^(k) 

/  n(u)du 
0 


y.(k) 


r (u)n(u)du 


l.e. , 


■  •  I\ 

l  a  (a  -  a  )  J  a  (6 6.  ' 

1=1  1  1  1-1  i=l  1  1  '-1 


where  =  fl  =  0  .  Note  that 


1^1  ai(ai  “i-l)  ~  (al  "  a2)oti  +  (a2  "  a3)a2  +  - 


■  +  ak\  = 


i  _ 


where  A , 


ai-l  —  ®  f°r  i  =  i*2,  k-1  and  A, 


Hence 


B1  6k  k  Vi  k  A-Bi 

~1~  implies  [  —  2.  1  '~t' 1  *  which  proves  (1)  . 

i  k  1=1  ak  1=1 


L^J  || 
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4.  APPLICATION  OF  TOTAL  TIME  ON  TEST 

Assuming  an  exponential  distribution,  the  results  of  Bickel  and  Doksum  (1968) 
may  be  used  to  establish  the  asymptotic  normality  of  W  in  the  incomplete  data 
case  for  selected  vectors  a  -  (a^,  a^)  .  Perhaps,  the  most  useful  test  is 

the  total  time  on  test  statistic.  In  the  case  of  a  complete  sample  of  size  n  , 

•k 

this  is  in  the  Bickel-Doksum  paper,  obtained  by  choosing  a^  ■  -i/(n  +  1)  , 

after  algebraic  manipulation.  In  the  case  of  incomplete  data  as  described  in  the 
introduction,  with  k  failures  observed;  the  total  time  on  test  statistic  is 


W  o 
a 


k-1  >(i) 

I  «■-«/ 

_ Z(l-» 


n(u)du 


i-1 


I  n(u)du 


> 


obtained  by  choosing  a°  ■  (k-l,k-2,  1,0)  . 

The  exact  distribution  conditioned  on  the  number  of  observed  failures  k  _>  2 
is  easily  computed  in  this  case.  Table  1  is  a  short  table  of  percentage  points. 
Note  that,  under 


W 

a 


O 


8-  IJX  + 


u2  + 


+  u 


k-1 


when  (i  «■  1,2,  k-1)  are  independent  uniform  random  variables  on  [0,1]  . 

k-1 

Since  the  distribution  of  Wa<>  is  symmetric  about  — —  ,  we  tabulate  upper 


percentiles  only. 
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Theorem  2. 


Let  $  be  a  monotone  test  of  versus  based  on  a  sample  of  Incomplete 

data  as  described  in  the  introduction.  Then 

/III)  ZM  \ 

E  $|  /  n(u)du,  ....  I  n(u)du|  |  F  IFR 

_  '  0  Z(k-1)  ' 


^  E[$(Y^,  Y^)  |  F  exponential)  , 


where  Y^,  Y^  independent  exponentially  distributed  random  variables. 


Proof : 

For  i  <  j  , 


zrU)  Ui) 

I  n (u)du  /  r(u)n(u)du 


.(  i  ) 

Z(i) 

n(u)du  / 

(i-1) 

Z(i-1) 

Yi 


r  (u)n(u)du 


The  inequality  follows  from  Lemma  3;  the  stochastic  equality  follows  from  Lemma  1, 


Thus  ♦(Y1,  ....  Yk)  1 


/  III)  ilk)  \ 

}>|  I  n(u)du,  /  n(u)du|  .  The 

\  0  Z(k-1)  / 


conclusion  follows  by  taking  expectations. 
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